Astrophysics and Space Science 
DOI 10.1007/s»«»««-»»«-»«»«-» 



Westerlund 1: bound or unbound? 

Mengel, S • Tacconi-Garman, L.E. 



OO 

O 

o 

(N 



© Springer- Verlag ■ 



43 ■ 

Or 
I ' 

o ■ 



o 

oo 
o 



X 



Abstract 

The Galactic open cluster Westerlund 1 (Wd 1) rep- 
resents the ideal local template for extragalactic young 
massive star clusters, because it is currently the only 
nearby young cluster which reaches a mass of around 
~10 5 Mq. The proximity makes spatially resolved stud- 
ies of its stellar population feasible, and additionally 
permits direct comparison of its properties with mea- 
surements of velocity dispersion and dynamical mass 
for spatially unresolved extragalactic clusters. 

Recently, we published the dynamical mass estimate 
based on spectra of four red supergiants. We have now 
identified six additional stars which allow a determina- 
tion of radial velocity from the wavelength covered in 
our VLT/ISAAC near-infrared spectra (CO bandhead 
region near 2.29 /xm), improving statistics significantly. 
Using a combination of stepping and scanning the slit 
across the cluster centre, we covered an area which in- 
cluded the following suitable spectral types: four red 
supergiants, five yellow hypergiants, and one B-type 
emission line star. 

Our measured velocity dispersion is 9.2km/s. To- 
gether with the cluster size of 0.86 pc, derived from 
archival near-infrared SOFI-NTT images, this yields a 
dynamical mass of 1.5x10 5 Mq. Comparing this to the 
mass derived via photometry, there is no indication that 
the cluster is currently undergoing dissolution. 
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1 Introduction 

The concept that globular clust ers may still be formin g 
today dates back fifteen years (Ho ltzman et al.lll992l) . 
and ever since young massive star clusters (YMSCs) 
have been the focus of intense studies. First, it needed 
to be verified that they actually have the properties 
which will make them comparable to (part of) the glob- 
ular clusters seen today after a Hubble time of evolu- 
tion. 

Once it was clear that, at least regarding mass, size 
and concentration, they could indeed be globular clus- 
ters in the making ( Whitmore fc Schweizer 1995), sur- 
vival times became an issue. A fraction between 10% 
and around 90% of the clusters seems to dissolve in each 



decade of time (Fall et al. 2004: Mensrel et al. 


Bastian et al. 


2005: Whitmore et all 2007k At 



the high end of these destruction rates is in conflict 
with the assumption that clusters at an age of around 
10 Myr, which have typically survived tens of crossing 
times, are in virial equilibrium. 

However, only clusters which can be spatially re- 
solved into single stars can be used to verify their 
dynamical state directly. The well known and com- 
paratively numerous Galactic clusters with maximum 
masses of around 10 4 M Q are not suitable as templates, 
because at an age of around 10 Myr, they are not nec- 
essarily expected to be in virial equilibrium, owing to 
longer crossing times, and the increased impact of in- 
teraction with galactic environment. Therefore, Wd 
1 turned out to be the best-suited template in the 
Milky Way when it was discovered that it s mass is 
much higher than initially (Westerlund! 1961 ) thought, 



around xl0 5 M o (jClark et al.l 120051 ). Given its young 



age of around 5 Myr, there is a considerable probability 
to find it in the process of disruption. 

In an attempt to determine the dynamical state of 
Wdl, we conducted high resolution spectroscopy in the 
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near-infrared. That is the only wavelength range (be- 
cause of the high extinction suffered by the cluster) 
which simultaneously allows a high signal-to-noise ra- 
tio, high spectral resolution, and large spatial coverage. 



2 Observations and data reduction 



shortward of the CO bandhead. For the age determina- 
tion, we used a rough flux calibration as an additional 
constraint on the sp ectral types. Nothing has changed 
for those results since Mengel &: Tacconi-Garmani ( 2007) .| 
we refer to that publication for the details. 
The resulting spectra are shown in Fig. [21 



Observations were conducted on VLT-ANTU on the 
nights of March If and 12, 2006 and consisted of two 
parts: 1) scanning the slit across a region while in- 
tegrating, covering the whole field shown in Figure [TJ 
and 2) stepping in steps of 0'.'3, the slit width, across 
the smaller region indicated by the black box in Fig- 
ure [T] In both cases, the 2' long slit was oriented in IM- 
S-direction. The large spatial coverage of the scan was 
achieved by performing scans of three different slit po- 
sitions, offset in N-S-direction by 108" each time, which 
allows for some spatial overlap. 

The central wavelength was set to be 2.31 /zm, and 
the actual wavelength range covered is 2.249 - 2.373 /im 
at a resolution of R~9000 (the maximum resolution 
ISAAC allows). 

Reduction of the data was performed using the IRAF 
packagS It included sky subtraction (a scan across an 
empty sky field 10' away was performed with identical 
integration time of 300 s directly after the cluster scan, 
and sky positions were interleaved in the stepping ob- 
servation every 10 minutes), bad pixel correction, and 
flatfielding. 

Wavelength calibration combined the information 
of the XeAr lamp sp ectra and the night sky lines 



_SP' 



( Rousselot et al.l 12000) . Spatial distortion was easily 



corrected due to the many star spectra along the slit. 

We extracted the individual spectra by averaging all 
the individual frames which included a star of inter- 
est, and then extracting from the combined frame by 
defining the aperture by hand. 

For the correction of the telluric absorption features, 
we used the hot, featureless stars covered by the stepped 
spectra (for example, W42a and a star without a known 
designation) . This gave a much better result than using 
the dedicated solar type telluric standards taken after 
the scans, which left some residuals around the Mgl 
line. 

Finally, we normalized all spectra which were used 
for the analysis by dividing by the continuum level de- 
termined from a small window (roughly 50 pixels) just 



1 IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



3 Analysis and Results 

The four brightest stars are assigned spectral types 
ranging between M3I and M6I (foreground objects are 
ruled out also because of the small spread in radial ve- 
locity) . 

We determined the rms radial velocity dispersion of 
all stars where this was possible from the wavelength 
range covered in our observations, which limited the 
selection to four red supergiants, five (assumed) yellow 
hypergiants (YHGs), and one sgB[e] type emission line 
star. See Table [T] for the individual radial velocities, 
and also some other properties like designation, spectral 
type, etc. 

A challenge, compared to the determination of ve- 
locity dispersion from just the RSGs, was the fact that 
the YHGs and the RSGs have almost no absorption 
lines in common at this wavelength range, and we also 
(currently) have no access to theoretical high-resolution 
spectra in the K-band for the YHGs. Therefore, we 
used a multi-step-procedure in an attempt to use a 
common wavelength reference frame for all our stars 
(except W9, which is a special case): First, we cross- 
correlated the spectrum of W237 in the region of CO 
absor ption with th e rest w avelength Arcturus spectrum 



from Hinklc et al 



( 19951 ). Cross-correlation of W237 
with the other RGSs around the CO absorption region 
gave us their radial velocities relative to W237. Using 
W20 and W26 as reference stars, we used the Ca absop- 
tion lines between 2.26 and 2.263/xm, and the Mg line 
around 2.282/Ltm to determine the relative velocites of 
those stars with respect to W32. And finally, we used 
five different lines/regions (all listed in table [1]) for the 
determination of relative velocities between W32 and 
the other YHGs. The radial velocity of W9 was deter- 
mined by manual best fit of the S emission lines. 

The resulting velocity dispersion is 9.2±2.5km/s. 
Uncertainties of this value were determined through 
Monte Carlo simulations, using the uncertainties in 
the individual velocity determinations, but we added 
0.5km/s uncertainty because of the additional uncer- 
tainty in matching RSG and YHG velocity reference 
points. 

This value is somewhat higher (even though consis- 
tent within the uncertainties) than the value we had ob- 



tained previously from the RSGs alone (jMengel fc Tacconi-Garma: 



Westerlund 1: bound or unbound? 



Star Reference star Spectral lines used v corr Spectral Type 



W26 


W237 


CO 


-48.7 


M5-6I 


W20 


W237 


CO 


-49.2 


M5I 


W75 


W237 


CO 


-59.2 


M4I 


W237 


Arcturus 


CO 


-61.7 


M3I 


W4 


W32 


lines below 2.26^m, S 2.27, Mg 2.28, C2.29 


-57.8 


F2Ia+ 


W8a 


W32 


lines below 2.26/rai, S 2.27, Mg 2.28, C2.29 


-41.9 


F5Ia+ 


W9 


W32 


S emission 


-62.4 


sgB[e] 


W12a 


W32 


lines below 2.26/im, S 2.27 


-61.3 


A5Ia+ 


W265 


W32 


lines below 2.26^m, S 2.27, C2.29 


-35.5 


F5Ia+ 


W32 


W20, W26 


Ca 2.26, Mg 2.28 


-57.4 


F5Ia+ 



Table 1 Selected properties of the stars used for the measurement of the velocity dispersion: Designation as in Wes terlundl 
l| 19871 ). the star which was used as a crosscorrelation reference, the spectral lines used for the cross correlation, the he- 
liocentrically corrected radial velo city (we applied a correcti on of +27.6km/s), and the spectral types (for the RSGs from 
iMengel fc Tacconi-Garmanl (|2007t ) , for the other stars from IClark et al.l (|2005h . See text for an explanation of reference 
stars and spectral lines. 
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Fig. 1 NTT/SOFI near- infrared image of Wd 1, composed 
from narrow band images (2.09/im and 2.17/im). The area 
shown was covered by the slit scans (with the slit oriented 
N-S), and the area shown by the black box was covered by 
steps of Of! 3 separation. Different stellar types are indicated 
by di fferent coloured cir cles, and the stars are identified by 
their IWesterlundl l|l987t ) designation. 



2007). This could have different reasons: The most 
obvious one would be simply statistics. Using only 
four stars before, statistic uncertainties are severe. The 
new value should be somewhat more reliable in that 
respect, but here the nature of the additional stars 
could i ntroduce a dditio nal velocity dispersion. Accord- 
ing to de Jager (1998), YHGs undergo radial pulsa- 
tions, which can cause variations in radial velocity of 
up to rs ±7km/s. Literature on spectral variations of 
YHGs is scarce, but from what is known, it seems easily 
possible that an increase in velocity dispersion of a few 
km/s could originate from this. 

The other parameter we require for the determi- 
nation of the dynamical mass is the projected half- 
light radius. We used the NIR image we created 
from the archival NTT/SOFI data to obtain a radius 
which contains half of the total light (in projection) 
to be r/ lp =0'.'86±0'.'14. The uncertainty comes mostly 
through assumption of different locations of the weakly 
constrained cluster centre. This value is consistent with 
the half-mass radius of lpc ( for stellar masses betwe en 
3.4 and 32 M Q ) presented bv lBrandner et all (|2007l ). 

The resulting dynamical mass, for a cluster in Virial 
equilibrium , is (assuming n = 9.3, for details on the 
formula see iMengel et ail(|2002f i) 



M dyn = W 2 r hp /G = 1.5lg;? x !0 5 M Q 

Within the uncertainties, this is consistent with that 
found bv lClark et all rt20Q5h. and somewhat high er than 
the mass determined by Brandner et al.l ( 2007t ). Both 
use stellar counts of the current population (with detec- 
tion of the main sequence only in the latter publication) 
to infer an initial cluster mass using an assumed stellar 
initial mass function. 
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Fig. 2 Specta of integrated cluster and the near-infrared brigthest individual stars. Atmospheric absorption (corrected 
for, but over-compensation and sub-op timal sky sub t ractio n causes a few artefacts at long wavelengths) is shown at the 
bottom. Line identifications taken from lHinkle et all l|l995t ). 



In iMengel fe Tacconi-Garmanl (|2007l ). we used the 
relative numbers of the stellar popul ations, together 
with evolutionary synthesis mo dels, (ILeitherer et al 



19991 Vazquez fe Leitherer 20051 ) to age date the clus- 



ter. We determine the best fitting age to be 5 Myrs. 

We used different methods (with reasonable agree- 
ment) to determine Ly for the cluster: Using the in- 
tegrated, extinction corrected K-band magnitude and 
the Starburst99 V-K colour at 5 Myr, or the integrated 
V-band magnitude with (substantial) extinction correc- 
tion, or the luminosity expected from Starburst99 for 
the cluster mass of 10 5 M Q expected from stellar counts. 
A total luminosity of Ly w 5.8 x 10 6 L Q leads to a light- 
to-mass-ratio Of Ly/M^yn ~ 40. 

Having ad ded our data point in the dia gnostic Figure 
[3j taken from Goodwin fc Bastianl < 20061 ) . we find that 
Wd 1 is located between the model lines denoting a star 
formation efficiency of 40% and 5 0%. Both of these ef- 
ficienc ies, in the N-body models of lGoodwin fc Bastian 
( 2006f h lead to the formation of stellar clusters which 
eventually, after 10-15 Myr, re-establish virial equilib- 
rium. This means that, based our data, Wd 1 does not 
seem to be currently in the process of dissolution. The 
potential influence of radial pulsations in the YHG sug- 
gests that the velocity dispersion and hence the virial 
cluster mass may be overestimated. On the other hand, 



mass segregation could have biased our measurements, 
both the velocity dispersion, and the cluster size, to- 
wards values which are too low. 



4 Conclusions and future work 

1. We identify ten stars which allow a determination of 
radial velocity in the part of the K-band which we 
use for our spectral analysis: four RSGs, five YHGs, 
and one emission line star. The rms velocity dis- 
persion of these stars was found to be 9.2±2.5km/s. 
The resulting dynamical mass is 1.5±l% x 1O 5 M . 

2. There is no indication of cluster dissolution. Com- 
paring L/M to theoretical expectations suggests a 
star formation efficiency between 40 and 50%, which 
is compatible with the cluster re-establishing equilib- 
rium after the phase of gas expulsion. 

3. We cannot rule out systematic or statistical errors 
in our analysis yet, since the sample of stars, even 
though larger than before, with just ten stars is still 
quite small. Furthermore, our measurements may 
have suffered a bias from either pulsations of the 
YHGs (velocity dispersion, and hence mass, would 
be too high) , or by not taking into account potential 
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mass segregation (velocity dispersion and size, and 
hence mass, would be too low). 

It is rather straightforward to overcome the poten- 
tial problems mentioned: For better statistics, we need 
to have access to more radial velocity determinations. 
We were granted observing time to observe other stel- 
lar types in a different near infrared wavelength region. 
Furthermore, we will observe two additional epochs, 
which will give us an indication of the impact of radial 
pulsation in the YHGs. And finally, there are two w ays 
to address mass segregation: iBrandner et al.l (|2007l ) see 
some indication of mass segregation in their SOFI data. 
Once deeper NIR imaging data, which already exist, 
are analyzed, an even better estimate of this property 
will be possible. Our new spectra will also cover lower 
mass stars, which will allow us to check if their velocity 
dispersion is lower than that of the lower mass stars. 

All of these future investigation will help settle 
the question conclusively if Westerlund 1 will dissolve 
within the next 10-20 Myr or not - an important step to- 
wards understanding what fraction of the galactic stel- 
lar population originates in star clusters. 



5 Questions and answers 

Mark Gieles: Perhaps there is some evidence for mass 
segregation: If the high mass stars sink to the centre, 
they will speed up, because they move deeper in the po- 
tential. Of course, this argument works only if you can 
confirm that your massive stars are also more centrally 
concentrated. 

Sabine Mengel: At the time of the presentation, deep 
NIR imaging data had not been published yet, and 
therefore the main sequence had been undetected and 
no information regarding mass segregation was avail- 
able. In the meantime, IBrandner et al.l (|2007l ) pub- 
lished their SOFI images, and they believe that they see 
some indication for mass segregation, however it is not 
quite clear for me yet if this conclusion is unaffected by 
statistical and/or observational bias. It is however be- 
lieved that the high mass stars segregate because they 
lose energy, and that they move more slowly than the 
low mass stars, once they are segregated. 
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Fig. 3 Evolution of L/Md yn for clusters with variable star 
form ation efficiencies (mo d ified version of a diagram taken 
from iGoodwin fc Bastiar] (|2006h V Clusters with an effec- 
tive star formation efficiency of below ss40% dissolve after 
around 10 Myr. Wdl (the black cross) is located in a re- 
gion which would correspond to a star formation efficiency 
of around 45%, and hence appears stable against (rapid) 
dissolution. 
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